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Abstract. We prove that the sequence (e®)„>o of excursions in the quarter plane corresponding 
to a nonsingular step set (5 C {0, ±1} 2 with infinite group does not satisfy any nontrivial linear re- 
currence with polynomial coefficients. Accordingly, in those cases, the trivariate generating function 
of the numbers of walks with given length and prescribed ending point is not D-finite. Moreover, 
we display the asymptotics of e®. 



1. Introduction 

1.1. General context. Counting walks in a fixed region of the lattice Z rf is a classical problem in 
probability theory and in enumerative combinatorics [48, 52, 26, 47, 46, 19, 16, 27, 56, 57, 6, 38, 5]. 
In recent years, the case of walks restricted to the quarter plane N 2 = £ Z 2 | i > 0, j > 0} 
has received special attention, and much progress has been done on this topic [35, 41, 32, 31, 24, 13, 
12, 11, 43, 9, 34, 44, 45, 14, 22, 18, 23, 37, 36, 25, 49, 7]. Given a set 6 of allowed steps, the general 
problem is to study ©-walks in the quarter plane N 2 , that is walks confined to N 2 , starting at (0, 0) 
and using steps in © only. Denoting by fe(i,j,n) the number of such walks that end at and 
use exactly n steps, the main high-level objective is to understand the generating function 

F e (x,y,t)= fe{i,3,n)xYt n e®[[x,y,t}}, 

i,j,n>0 

since this continuous object captures a great amount of interesting combinatorial information about 
the discrete object fe(i,j,n). For instance, the specialization Fe(l,l,t) counts ©-walks with 
prescribed length, the specialization i*e(l,0,i) counts ©-walks ending on the horizontal axis, and 
the specialization Fe(0, 0, i) counts ©-walks returning to the origin, called ©- excursions. 

1.2. Questions. From the combinatorial point of view, the ideal goal would be to find a closed 
form expression for /©(«,j, n), or at least for F<s(x,y,t). This is not possible in general, even 
if one restricts to particular step sets ©. Therefore, it is customary to address more modest, 
still challenging, questions such as: determine the asymptotic behavior of the sequence f&{i,j,n)\ 
determine the structural properties of F&(x,y,t): is it algebraic (that is, root of a polynomial in 
Q[x,y,t,T])7 is it D-finite? (in one variable t this means solution of a linear differential equation 
with coefficients in Q[i]; in several variables the appropriate generalization [40] is that the set of all 
partial derivatives spans a finite-dimensional vector space over Q(x,y,t)). These two questions are 
related, since the asymptotic behavior of the coefficient sequence of a power series is well understood 
for algebraic and D-finite power series [29]. 

1.3. Main result. In this work, we prove that the generating function 

F 6 (0,0,t) = J2 e nt n GQI1 

n>0 

of the sequence (e®) n >o of ©-excursions is not D-finite for a large class of walks in the quarter 
plane. Precisely, this large class corresponds to all small step sets © C {0, ±1} 2 for which a certain 
group G(&) of birational transformations is infinite, with the exception of a few cases for which 
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Fq(0, 0, t) = 1 is trivially D-finite (these exceptional cases are called singular, see below.) If % — X& 
denotes the characteristic polynomial of the step set S defined by 

(i,i)ee 

then the group G((3) is defined [24, 14] as a group of rational automorphisms of Q(x,y) that leave 
invariant the (Laurent) polynomial x( x :U)- Up to some equivalence relations, there are 51 cases of 
nonsingular step sets in N 2 with an infinite group. They are depicted in Table 1 in Appendix 3. 
With these definitions, our main result can be stated as follows. 

Theorem 1. Let & C {0, ±1} 2 be any of the 51 nonsingular step sets in N 2 whose group is infinite. 
Then the generating function _F@(0, 0, i) of &- excursions is not D-finite. Equivalently, the sequence 
of excursions (e®) n >o does not satisfy any nontrivial linear recurrence with polynomial coefficients. 
In particular, the full generating function Fe(x,y,t) is not D-finite. 

1.4. Previous results. 

1.4.1. Structural properties. While it is known that unrestricted walks in Z 2 have rational gener- 
ating functions and walks restricted to a half-plane in Z 2 have algebraic generating functions [3], a 
first intriguing result about walks in the quarter plane is that their generating functions need not 
be algebraic, and not even D-finite. For instance, Bousquet-Melou and Petkovsek [12] proved that 
this is the case for the so-called knight walk with 6 = {(2, —1), (—1, 2)}. (This actually constitutes 
one of the initial motivations to the study of walks evolving in the quarter plane.) It was later 
shown [45] that this remains true even if one restricts to next nearest neighbor walks: there still exist 
step sets © C {0,±1} 2 , for instance & = {(—1, 1), (1, 1), (1, —1)}, such that the series F@(x,y,t) is 
not D-finite. 

We restrict in the remaining of this text to next nearest neighbor walks, and reserve the wording 
small step walk to such particular walks in the quarter plane. Several sporadic cases of small 
step walks with D-finite generating functions have been known for a long time; among them, 
the most popular ones are Kreweras' walks [35, 11] for which 6 = {(— 1, 0), (0, — 1), (1, 1)}, and 
Gouyou-Beauchamps' walks [32] for which (3 = {(1, 0), (— 1, 0), (— 1, 1), (1, — 1)}. A whole class 
of small step walks with D-finite generating functions was first identified in [12, §2]: this class 
contains step sets that admit a vertical (or a horizontal) symmetry. Another class, including 
S = {(0,1), (-1,0), (1,-1)} and S = {(0,1), (-1,0), (1,-1), (0,-1), (1,0), (-1,1)}, corresponds 
to step sets that are left invariant by a Weyl group and whose walks are confined to a corresponding 
Weyl chamber [31]. 

A first systematic classification of small step walks with respect to D-finiteness was then un- 
dertaken by Mishna [43, 44], but only for step sets of cardinality at most three. A complete, still 
conjectural, classification was obtained in [9] using computer algebra tools. Almost simultaneously, 
Bousquet-Melou and Mishna [14] rigorously proved that among the 2 8 possible cases of small step 
walks there are exactly 79 inherently different cases of walks in the quarter plane, and they identified 
among them 22 cases of step sets & having a D-finite generating function FQ(x,y,t). 

A 23rd case, namely <8 = {(1,0), (—1,0), (1, 1), (—1, —1)}, known as Gessel walks, was discovered 
and proved to be D-finite, and even algebraic, in [10], using computer algebra techniques. It was 
proved afterwards in [22] by a different approach that for any fixed value to £ (0, 1/4), the bivariate 
generating function F<g{x,y,tQ) for Gessel walks is algebraic over R(x,y). To our knowledge, there 
is currently no "purely mathematical" proof of the algebraicity (and even of the D-finiteness) of 
the full generating function F&(x,y,t). 

Bousquet-Melou and Mishna [14] showed that the 23 cases of step sets & with D-finite generating 
function F&(x,y,t) correspond to walks possessing a finite group G(6). Informally speaking, the 
group of a walk is a notion that captures symmetries of the step set and that is used to generalize a 
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classical technique in lattice combinatorics called the "reflection principle" [26, Ch. III. 1] . Moreover, 
Bousquet-Melou and Mishna [14] conjectured that the 56 remaining models with infinite group have 
non-D-finite generating functions Fq(x, y, t). This was recently proved in [37] for the 51 nonsingular 
walks, that is, for walks having at least one step from the set {(— 1, 0), (— 1, — 1), (0, — 1)}. It is 
worth mentioning that for these 51 models, the non-D-finiteness of F<s(x,y,t) was obtained as a 
consequence of the non-D-finiteness of this series as a function of x, y. Further, before our work, 
nothing was known about the non-D-finiteness of -Fg(0,0, i). Two out of the five singular walks 
were already shown to have non-D-finite generating series in [45]. Therefore, before this work, there 
still remained 3 cases of singular walks for which the generating function Fq(x, y, t) was suspected, 
but not yet proved, to be non-D-finite [42]. 

1.4.2. Closed form expressions. Kreweras [35] and Gouyou-Beauchamps [32] found explicit counting 
formulas for the walks named after them. For Gessel walks, an explicit formula for excursions was 
conjectured around 2000 by Ira Gessel and proved by Kauers, Koutschan and Zeilberger in [34]. A 
closed form expression (in terms of nested radicals) for the full Gessel generating function F<$(x, y, t) 
was proved in [10]. Some other explicit formulas for fe(i,j,ri) and for F<s(x,y,t), or some of their 
specializations, have been obtained in cases when 6 admits a finite group [14]. A different type 
of explicit expressions (integral representations) for the generating function of Gessel walks was 
obtained in [36]. The approach of [36] was later generalized in [49] to all the 74 nonsingular walks. 
Finally, the article [7] provides unified formulas in terms of Gaussian hyper geometric functions for 
all D-finite transcendental series F(s(x,y,t). Again, no "purely human" (computer-free) proof of 
this result is known yet. 

1.4.3. Asymptotics. Concerning asymptotics, conjectural results were displayed in [9] for the coeffi- 
cients of Fq(1, 1, t) when this latter function is D-finite. Some of these conjectures have been proved 
in [14]. Explicit asymptotics for the coefficients of F&(0, 0,t) and F&(1, 1, i) were conjectured even 
in non-D-finite cases in some unpublished tables [8]. In a recent work, Denisov and Wachtel [18] 
obtained explicit expressions for the asymptotics of excursions Fq(0, 0, t) in a much broader setting; 
in particular, their results provide (up to a constant) the dominating term in the asymptotics of the 
n-th coefficient of Fq(0, 0,t) in terms of the step set. Even more recently, Fayolle and Raschel [25] 
showed that the dominant singularities of i*e(0, 0,i), i<e(l,0, t) and Fq(1, l,t) are algebraic num- 
bers, and announced more general and precise results about asymptotics of coefficients of Fq(0, 0, t), 
F 6 (l,0,i) andF 6 (l,l,i) [21]. 

2. Probability, Number Theory and Algorithms 

2.1. Contributions. In the present work, we prove the non-D-finiteness of the generating series 
of excursions Fq(0, 0, t) for the 51 cases of nonsingular walks with infinite group. As a corollary, 
we deduce the non-D-finiteness of the full generating function F<g(x,y,t) for those cases. Indeed, 
D-finite series are closed under Hadamard product [40] and 1/(1 — t) is clearly D-finite with respect 
to x,y,t. This corollary has been already obtained in [37], but the approach here is at the same 
time simpler, and delivers a more accurate information. This new proof only uses asymptotic 
information about the coefficients of F&(0, 0, t), and arithmetic information about the constrained 
behavior of the asymptotics of these coefficients when their generating function is D-finite. More 
precisely, we first give consequences of the general results in [18] in the case of walks in the quarter 
plane. If e n = e® denotes the number of excursions of length n using only steps in 6, this analysis 
implies that, when n tends to infinity, e n behaves like K ■ p n ■ n a , where K = K{&) > is a real 
number, p = p{&) is an algebraic number, and a = a(&) is a real number such that c = cos(j^) 
is an algebraic number. Moreover, explicit real approximations for p, a and c can be determined to 
arbitrary precision, and exact minimal polynomials of p and c can be determined algorithmically 
starting from the step set &. For the 51 cases of nonsingular walks with infinite group, this enables 

3 



us to prove that the constant a = a(6) is not a rational number. The proof amounts to checking 
that some explicit polynomials in Q[t] are not cyclotomic. To conclude, we use a classical result in 
the arithmetic theory of linear differential equations [20, 2, 30] about the asymptotic behavior of 
an integer-valued, exponentially bounded D-finite sequence, stating that if such a sequence grows 
like K ■ p n ■ n a , then a is necessarily a rational number. 

In summary, our approach brings together (consequences of) a strong probabilis- 
tic result [18] and a strong arithmetic result [20, 30], and demonstrates that this 
combination allows for the algorithmic certification of the non-D-finiteness of the 
generating series of excursions Fq(0, 0, t) in the 51 cases of nonsingular walks with 
infinite group. 

2.2. Number theory. It is classical that, in many cases, transcendence of a complex function can 
be recognized by simply looking at the local behavior around its singularities, or equivalently at 
the asymptotic behavior of its Taylor coefficients. This is a consequence of the Newton-Puiseux 
theorem and of transfer theorems based on Cauchy's integral formula, see, e.g., [28, §3] and [29, 
Ch. VII. 7]. For instance, if (a n ) n >o is a sequence whose asymptotic behaviour has the form K-p n -n a 
where either the growth constant p is transcendental, or the singular exponent a is irrational or 
negative integer, then the generating series A(t) = J2n>o a n^ n is n °t algebraic. 

A direct application of this criterion to our case allows to show at a glance that the generating 
series for excursions in the 51 cases of nonsingular walks with infinite group are transcendental. 

Similar (stronger, though less known) results, originating from the arithmetic theory of linear 
differential equations, also allow to detect non-D-finiteness of power series by using asymptotics of 
their coefficients. This is a consequence of the theory of G- functions [1, 20], introduced by Siegel 
almost a century ago in his work on diophantine approximations [51]. 

We will only use a corollary of this theory, which is well-suited to applications in combinatorics. 

Theorem 2. Let (a n ) n >o be an integer-valued sequence whose n-th term a n behaves asymptotically 
like K ■ p n ■ n a , for some real constant K > 0. If the growth constant p is transcendental, or if the 
singular exponent a is irrational, then the generating series A(t) = Xm>o a nt n is not D-finite. 

Proof. This result is more or less classical, but we could not find its exact statement in the literature. 

Classical results by Birkhoff-Trjitzinsky [4] and Turrittin [53] imply that if the n-th coefficient of 
a D-finite power series is asymptotic to K ■ p n ■ n a , then p and a are necessarily algebraic numbers. 

The difficult part of Theorem 2 is that irrationality of the singular exponent implies non-D- 
finiteness, under the integrality assumption of coefficients. The only proof that we are aware of 
uses the fact that any D-finite power series with integer-valued and exponentially bounded coef- 
ficients is a G-function. It relies on the combination of several strong arithmetic results. First, 
the Chudnovsky-Andre theorem [15, 1] states that the minimal order linear differential operator 
satisfied by a G-function is globally nilpotent. Next, Katz's theorem [33] shows that the global 
nilpotence of a differential operator implies that all of its singular points are regular singular points 
with rational exponents. 

We refer to [20] for more details on this topic, and to [30] for a brief and elementary account. □ 

2.3. Probability. Random processes in cones have arisen a great interest in the mathematical 
community, as they appear in several distinct domains: quantum random walks, random matri- 
ces, noncolliding random walks, etc. Accordingly, many results are available concerning the exit 
times from cones of the Brownian motion, and more recently for random walks in discrete time. 
For our purpose, the setting is the following: given a sequence (Xi(k), X2(k))k>i of independent 
and identically distributed (i.i.d.) random variables such that for all step s £ & and all k > 1, 
¥[(Xi(k), X2(k)) = s] = 1/|@ |, the exit time r is the hitting time of the boundary of the (translated 
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positive) quarter plane ({-1} UN) 2 , i.e., r = inf{n > 1 : J2k=i x i( k ) = - 1 or T% =1 X 2 (k) = -1}. 
Now, the standard relation between probability and counting reads 



(1) 



\6\ n 



J2(Xi(k),X 2 (k)) = (i,j),r>n 
Dc=l 

The asymptotic behavior of this probability, as the time n goes to infinity, is called a local limit 
theorem. This is now well understood for a very large class of random walks and cones [18, 55]. To 
the best of our knowledge, even in the special case of the small step nonsingular random walks in the 
quarter plane, the precise asymptotic behavior of (1) was not known before [18]. The corresponding 
result is the following. 

Theorem 3 ([18]). Let & C {0, ±1} 2 be the step set of a nonsingular walk in the quarter plane N 2 . 
Let e n = e® denote the number of excursions of length n using only steps in &, and let x = X& 
denote the characteristic polynomial J2(i,j)ee X% V^ ^ Q[ x i Ui V~ l ] °f the step set 6. Then 

- if the walk is aperiodic, 

e n ~K-p n - n a , 

- if the walk is periodic (then of period 2), 

e 2n ~K-p 2n -(2n) a , e 2n+1 = 0, 

where in both cases, K = Kq is a real constant, p = pe is an algebraic number, and a = a© is a 
real number such that cos(yq^) is an algebraic number. Moreover, the system 

^ ^ dx dy 

has a unique solution (xo,2/o) G ^>o- ^ ^ s suc -h that p = x{ x o-,Vo) an d a = —1 — it/ arccos(— c), 
where c is the algebraic number 

d 2 x 

(3) c= W=^ ( x ViV°)- 

V dx' 2 dy 2 

Note that actually, the precise lower and upper bounds for the probability in (1) obtained pre- 
viously by Varopoulos [55] (see in particular [55, Theorem 4]) would be sufficient for our needs, as 
they already provide the values of p and a. 

Proof of Theorem 3. This is a rewriting of results in [18] . Consider first an aperiodic step set © 
and a random walk with increments (Yi(k),Y2(k))k>i such that for all s € & and all k > 1, 

(4) F[(Y 1 (k),Y 2 (k)) = s] = eM J^ S) \ R(h) = $>xp«M), 

where h = (hi,h 2 ) is fixed so that E,[(Yi(k),Y 2 (k))] = 0. It turns out that for all nonsingular 
cases, this uniquely determines h (the uniqueness follows from the strict convexity of the Laplace 
transform R(h), see [57, Ch. 8]). With this notation, [18, Formula (9)] reads 

(5) fe(i,j,n) = C((0,0),(i,j))R(h) n n-P- d / 2 (l + O (l)), n -> oo, 
where 

- C((0, 0), is some constant (in fact, the product of certain harmonic functions evaluated 
at the starting and ending points (0,0) and 

- d is the dimension of the lattice where the random walk evolves, here equal to 2; 

- p = vr/arccos(-r), with r = E[y i (A;)y 2 (A;)]/ v /E[Yi(A;) 2 ] • E[Y 2 {k)' 2 }, see [18, Example 2]. 
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To prove Theorem 3, we thus have to show that R(h) = x( x o, 2/o)> and that r = c, with c as in (3). 

We start with proving that R(h) = x( x QiUo)- The law of the increments (Yi(k),Y 2 (k)) of the 
random walk has the expectation 

E[(Y 1 (k),Y 2 (k))} = 1 ±- £ (.i,j)eM((hi,h 2 ),(i,j)}). 
{ ' (i,i)ee 

Accordingly, E[(Y 1 (k),Y 2 (k))] = if and only if 

J2 i(eMhi)y(exp(h 2 )y = £ i(exp(/ l i)) i (exp(/ i2 )y = 0. 

This exactly means that (xo, yo) = (exp(hi), exp(/t 2 )) satisfies the system (2), and then that R(h) = 

xi x o,yo)- 

We now show that r = c, with c as in (3). Using (4) yields 

E\Yi(k)Y 2 (k)] = E(i, i )g 6 ^( ex P(^i)) t ( ex P(^)) i 

V«PM " v /(E( JJ ) e6 i 2 (exp(/ i i)) i (exp(/ i2 ))i) • (E( Jj)ee J 2 (exp(/ ll )) i (exp(/ l2 )y)' 

The numerator above is precisely xoyo-Q^(xQ,yo). As for the denominator, we notice that 

jr4(x ,y ) = J2 ^~ l ) x V 2 Vo = h *(* ~ 1 )4»o = i E * 2x oyo: 
(ij)ee x °(i,i)Ge x o (iji)e6 

where the last equality comes from ^(xo,yo) = 0. A similar computation for 2^(xo,yo) shows 
that r = c. 

We now consider the periodic step sets (they are marked with a star in Tables 1 and 2). By 
examining the different cases, it is clear that the period is necessarily 2. In particular, we have 
C2n+i = for all n > 0, and (5) does not hold anymore. However, a slight adaptation of the local 
limit theorem [18, Theorem 6] (from which the asymptotic formula (5) is a simple consequence) 
gives 

fe(i,j,2n) = C((0,0), (i, j))R(h) 2n (2 n y^ d / 2 (l + o(l)), n -> oo, 

where C((0, 0), (i, j)) is some constant (different from that in (5)). The proof of Theorem 3 is 
complete. □ 

2.4. Algorithmic irrationality proof. Let (3 C {0, ±1} 2 be one of the 51 nonsingular step sets 
with infinite group (see Table 1 in Appendix 3). By Theorem 3, the singular exponent a in the 
asymptotic expansion of the excursion sequence (e®) n >o is equal to —1 — it/ arccos(— c), where c 
is an algebraic number. Therefore, if arccos(c)/-7r is an irrational number, then by Theorem 2, the 
generating series F&(0,0,t) is not D-fmite. 

We now explain how, starting from the step set 6 one can algorithmically prove that arccos(c)/7r 
is irrational. This effective proof decomposes into two main steps, solved by two different algorithms. 
The first algorithm computes the minimal polynomial fA c (t) € Q[t] of c starting from &. The second 
one performs computations on // c (i) showing that arccos(c)/7r is irrational. 

2.4.1. Computing the minimal polynomial of the correlation coefficient. Given x = Xe the char- 
acteristic polynomial of the step set 6, Theorem 3 shows that the exponential growth p and the 
correlation coefficient c are algebraic numbers, for which algebraic equations can be obtained by 
eliminating x and y from the equations 



This elimination is a routine task in effective algebraic geometry, usually performed with Grobner 
bases for lexicographic or elimination orders [17]. For any zero (xo,yo) of the system §f (xo,yo) = 
^(xQ,yo) = and any polynomials P(x,y) and Q(x,y) such that Q $ I, the algebraic num- 
ber P(xq, yo)/Q(xo,yo) is a root of a generator of the ideal / + (P(x, y) — tQ(x, y)) n Q[t\. 
This is summarized in the following algorithm. 

Input: A step set & satisfying the assumptions of Theorem 3 
Output: The minimal polynomials of p and c defined in Theorem 3 

(1) Set x(x,y) ■■= E(i,j)ee S V > and compute Xx ■= numer(^), x y ■= numer(§|). 

(2) Compute the Grobner basis of the ideal generated by (xx, Xy, numer(t — x)) m Q[x, y, t] for 
a term order that eliminates x and y. Isolate the unique polynomial in this basis that is 
free of x and y, factorize it, and identify its factor p p that annihilates p. 

(3) Compute the polynomial 



and eliminate x and y by computing a Grobner basis of {XxiXyiP) f° r a term order that 
eliminates x and y. Isolate the unique polynomial in this basis that is free of x and y, 
factorize it, and identify its factor p c that annihilates c. 

Table 2 in Appendix 3 displays the minimal polynomials of p and of c obtained using this 



2.4.2. Proving that the arccosine of the correlation coefficient is not commensurable with tt. Given 
the minimal polynomial p c of the correlation coefficient c, we now want to check that arccos(c)/7r is 
irrational. General classification results exist, e.g., [54], but they are not sufficient for our purpose. 
Instead, we rather prove that arccos(c)/7r is irrational in an algorithmic way. This is based on the 
observation that if arccos(c)/7r were rational, then c would be of the form (x + l/x)/2 with x a 
root of unity. This implies that the numerator of the rational function fj, c ( x J~ ) would possess a 
root which is a root of unity. In other words, the polynomial R{x) = x de z^fj, c (^) would be 
divisible by a cyclotomic polynomial. This possibility can be discarded by analyzing the minimal 
polynomials p c displayed in Table 2 in Appendix 3. 

Indeed, in all the 51 cases, the polynomial R(x) is irreducible and has degree 2deg(/i c ), thus at 
most 28. Now, it is known that if the cyclotomic polynomial $at has degree at most 30, then is at 
most 150 [50, Theorem 15], and the coefficients of &n belong to the set {—2, —1, 0, 1, 2} [39]. Com- 
puting R in the 51 cases shows that it has at least one coefficient of absolute value greater than 3. 
This allows to conclude that R is not a cyclotomic polynomial, and therefore that arccos(c)/7r is 
irrational, and finishes the proof of Theorem 1. 

2.5. Example. We now illustrate the systematic nature of our algorithms on Example 23 of Ta- 
ble 1, i.e., walks with step set 6 = {(—1, 0), (0, 1), (1, 0), (1, —1), (0, —1)}. For ease of use, we give 
explicit Maple instructions. 

Step 1. The characteristic polynomial of the step set is 

S: = [[-1,0] , [0,1] , [1,0] , [1,-1] , [0,-1]] : 
chi:=add(x~s[l]*y~s[2] ,s=S) ; 




algorithm. 
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11 X 

X-= -H hx + y+-, 

x y y 

whose derivatives have numerators 

chi_x: =numer (dif f (chi ,x) ) ; chi_y :=numer (dif f (chi ,y) ) ; 

(6) Xx ■= x 2 + x 2 y - y, x y ■= y 2 ~ x - 1. 

These define the system (2). 

Step 2. We now compute a polynomial that vanishes at p = x( x o,yo) when (#o,yo) is a solution 
of (6). To this aim, we eliminate x and y in {x x , Xyi numer(x) — tdenom(x)} by a Grobner basis 
computation using an elimination order with (x, y) > t. In Maple, this is provided by the command 

G : =Groebner [Basis] ( [chi_x, chi_y ,numer (t-chi)] ,lexdeg( [x,y] , [t] )) : 

which returns four polynomials, only one of which is free of x and y, namely 

p : =f actor (op (remove (has , G , {x , y}) ) ) ; 

p:= (t+l)(t 3 + t 2 -18t-43). 

Since we know that p > 0, we identify its minimal polynomial as p p = t 3 + t 2 - 18t - 43, which 
gives the entry in Column 3 of Table 2. The numerical value for p in Table 1 is given by 

f solve (p,t ,0 . . infinity) ; 

4.729031538. 

Step 3. Next, we obtain a polynomial which vanishes at c by a very similar computation: 

G : =Groebner [Basis] ( [chi_x , chi_y , 

numer (t~2-dif f (chi ,x,y) ~2/dif f (chi ,x,x) /dif f (chi ,y ,y) )] , lexdeg( [x,y] , [t] ) ) ; 

Again, this command returns four polynomials, with one of them free of x and y, namely 

p : =f actor (op (remove (has , G , {x , y}) ) ) ; 

p := {At 2 + l)(8i 3 + 8i 2 + 6t + l)(8i 3 - 8t 2 + 6t - 1). 

This polynomial has only two real roots, ±c. Since c < 0, we identify its minimal polynomial 
as p c = 8t 3 + 8t 2 + 6t + 1, which gives the entry in Column 4 of Table 2. Again, the numerical 
value for a in Table 1 is given by 

evalf (-1-Pi/arccos (-f solve (p,t , -infinity . . 0) ) ) ; 

-3.320191962. 

Step 4- To conclude, we compute the polynomial 

(x 2 + 1 \ 
— J = x 6 + 2x 5 + 6x 4 + 5x 3 + 6x 2 + 2x + l. 

This polynomial does not have any root that is a root of unity, since it is irreducible and not 
cyclotomic: 

irreduc(R) ,numtheory [iscyclotomic] (R,x) ; 

true, false 

This completes the proof that the generating function for this walk is not D-finite. 

2.6. Open problems. As we mentioned in Section 2.1, our approach brings together a strong 
arithmetic result (Theorem 2) and a strong probabilistic result (Theorem 3). It therefore appears 
natural to search for alternative simpler proofs of these results. 



Proving that a is transcendental. In Section 2.4, we are able to prove that for the 51 nonsingular 
models, the exponent a in the asymptotic expansion of the excursion sequence is irrational. It is 
worth mentioning that if it was possible to prove that a is not only irrational, but also transcen- 
dental, then Theorem 2 would not be needed. 

Simpler proof of Theorem 2. In Section 2.2, we proposed a proof of Theorem 2 based on several 
strong results from arithmetic theory [15, 1, 33, 20]. It would be interesting to know whether 
Theorem 2 admits a simpler, direct proof. 

Combinatorial proof of Theorem 3. Last but not least, finding a combinatorial proof of Theorem 3 
appears as an interesting challenge. 
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3. Appendix 



Table 1 Non-singular walks in the quarter plane with an infinite { 
their excursions. The labels used in the columns "Tag" correspond 
Entries marked with a star correspond to periodic walks. 



;roup, and the asymptotics of 
to those from Table 4 in [14]. 



Tag Steps 



First terms 



Asymptotics Tag Steps 



First terms 



5 
7* 
9 

11* 
14 
17* 
19* 
21 
24 
26 
28 
30 
32* 
34* 
36 
38 
40 
42 
44 
46 



7^ 



1, 0, 1, 2, 2, 13, 21, 67, 231 



1, 0, 1, 0, 4, 0, 29, 0, 230 



1, 0, 1, 1, 2, 7, 10, 38, 89 



1, 0, 0, 0, 2, 0, 6, 0, 42 



1, 0, 0, 1, 2, 0, 5, 26, 28 



1, 0, 1, 0, 4, 0, 29, 0, 230 



1, 0, 0, 0, 2, 0, 6, 0, 42 



1, 0, 1, 1, 4, 7, 25, 64, 201 



1, 0, 2, 2, 10, 26, 86, 312, 1022 



1, 0, 2, 2, 11, 27, 101, 348, 1237 



1, 0, 1, 2, 4, 13, 36, 111, 343 



1, 0, 1, 1, 6, 17, 58, 202, 749 



1, 0, 2, 0, 13, 0, 124, 0, 1427 



1, 0, 1, 0, 5, 0, 35, 0, 313 



1, 0, 0, 1, 2, 1, 11, 27, 60 



1, 0, 2, 2, 10, 26, 86, 312, 1022 



1, 0, 0, 2, 4, 8, 28, 108, 372 



1, 0, 0, 2, 4, 8, 28, 108, 372 



1, 0, 2, 3, 15, 51, 208, 893, 3841 



1, 0, 1, 2, 8, 22, 101, 364, 1618 



3.799605" 

n 2. 610604 

3.799605™ 



n 



2.318862 



3.800378™ 



3.799605™ 



3.800378" 

n 3. 918957 

3.608079" 

n 3. 388025 

3.800378" 

n 2. 521116 

3.800378" 

n 3. 918957 

4.214757" 



4.729032" 

n 2. 757466 

4.729032" 



2.397625 



4.214757" 



4.569086" 

n 2. 503534 

4.569086" 

n 3. 985964 

4.372923" 



4.070925 



4.729032" 

^2. 757466 



5.891838" 



5.799605" 



4 

6 

8 
10 
12 
10 
18 
20 
23 
25 
27* 
29* 
31 
33 
35 
37 
39 
41 
43 
45 
47 



1, 0, 0, 2, 2, 0, 16, 44, 28 



1, 0, 1, 2, 2, 13, 21, 67, 231 



1, 0, 1, 1, 2, 7, 10, 38, 



1, 0, 0, 1, 2, 0, 5, 26, 28 



1, 0, 0, 1, 0, 1, 5, 1, 18 



1, 0, 1, 2, 2, 14, 21, 76, 252 



1, 0, 0, 2, 2, 0, 16, 44, 28 



1, 0, 1, 2, 4, 14, 45, 120, 468 



1, 0, 2, 1, 10, 14, 75, 178, 738 



1, 0, 2, 2, 11, 27, 101, 348, 1237 



1, 0, 2, 0, 13, 0, 124, 0, 1427 



1, 0, 1, 0, 5, 0, 35, 0, 313 



1, 0, 0, 1, 2, 1, 11, 27, 60 



1, 0, 1, 1, 4, 7, 25, 64, 201 



1, 0, 1, 1, 3, 8, 19, 65, 177 



1, 0, 1, 2, 4, 13, 36, 111, 343 



1, 0, 1, 1, 3, 8, 19, 65, 177 



1, 0, 1, 2, 4, 14, 45, 120, 468 



1, 0, 2, 2, 13, 27, 140, 392, 1882 



1, 0, 1, 1, 5, 8, 40, 91, 406 



1, 0, 1, 3, 7, 29, 101, 404, 1657 
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Tag Steps 



First terms 



Asymptotics Tag Steps 



First terms 



Asymptotics 



50 
52 
54 
56 



1, 0, 1, 1, 5, 8, 40, 91, 406 
1, 0, 2, 3, 15, 51, 208, 893, 3841 
1, 0, 1, 1, 8, 18, 90, 301, 1413 
1, 0, 3, 5, 30, 111, 548, 2586, 13087 



5.064419" 

n 4. 036441 

5.891838" 

n 2. 679783 

5.799605" 



6.729032" 



6.729032" 

„3. 497037 



49 
51 

53 
55 



1, 0, 2, 2, 13, 27, 140, 392, 1882 
1, 0, 1, 3, 7, 29, 101, 404, 1657 
1, 0, 1, 2, 8, 22, 101, 364, 1618 

1, 0, 2, 4, 16, 64, 266, 1210, 5630 



5.064419" 

n 2. 491053 

5.891838" 

n 3. 471058 

5.799605" 



6.729032" 



Table 2 Non-singular walks in the quarter plane with an infinite 
mials of the growth constants p and of the correlation coefficients 



Stops 



Minimal polynomial /j, p of p 



group, and the minimal polyno- 
c. 

Minimal polynomial fi c of c = — cos( _7 ) 



f 4 + t 3 - 8i 2 - 36t - 11 



,8 , 1 ,6 _ 3 ,4 , 3 +2 _ 1 



1 ,6 _ 3 .4 , 3 ,2 _ 1 



,4 _ 9,3 , 27 ,2 _ 35 , , 17 
1 2 L " p 4 1 8 ^ 16 



*4 9,3 , 27 ,2 I 35 f , 17 



t 6 - lit 4 - 32t 2 - 256 



t 6 + | t 4 + 2t 2 _ 1 



t 6 + ft 4 + 2t 2 - | 



t 5 + f 4 + i 3 - 30i 2 - 96t - 91 



! + t 6 - 


64 1 


- -^i 2 - 

256 


1 

1024 


! + t 6 - 


it 4 - 
64 1 


- -M 2 - 

256 


1 

1024 


157 ,6 , 
32 1 "+" 


145 ,4 
128 


, 1681 .2 
" + " 512 1 


2209 
2048 


157,6 , 
32 1 "I" 


145 ,4 


, 1681 ,2 
"t" 512 1 


2209 
2048 



i 5 - 2t 4 - At 3 - 31t 2 + 23t - 41 



i 5 + 2t 4 - 7i 3 - 46i 2 - 116t - 131 



,10 , 3 ,8 , 13,6 , _5_,4 , _3_,2 1_ 

L 2 1 16 "i" 64 1 ^ 256 1024 



,10 _i_ 3,8 _i_ 13,6 _i_ _5_/4 j 3_,2 _ 1 

_ 2 1 - ^ IC 1 l_SAi: 1 256 1024 



t 3 + t 2 - 18t - 43 



t Z + t 2 + ft + i 



t 3 _ j2 + | t „ 1 



+4 , 7,2 _ 5 



f6 _ fi r JL t 2 _ _5_ 



t 6 - 20t 4 - 16i 2 - 48 



t 6 + 2t 4 + 1 t 2 - | 
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Tag 


Steps 


Minimal polynomial /i p of p 


Minimal polynomial fi c of c = — cos( ) 


30 
40, 42 


WW 


t - 5 


*+* 


43, 49 
45, 48 


W\ 




t 6 + 2t 5 - 18i 4 - 67t 3 - 108t 2 - 40i - 19 


,12 i 11,10 , 107,8 , 145,6 , 455,4 2859,2 , 1521 
~ r 4 ' 16 ' 32 ' 128 1024 ' 4096 

,12 , 11,10 , 107,8 , 145,6 , 455,4 2859,2 , 1521 
4 1 16 32 1 128 1 " 1024 1 - 1" 4096 


44, 50 
47, 51 


W , j/ 

WW 


t 7 + 3t 6 - 18t 5 - 127t 4 - 328t 3 - 560i 2 - 704t - 448 


,14 , 23 ,12 , 25,10 i 971,8 i 421 ,6 , 307,4 , 107,2 49 
'4 ~ r 2 ~ r 64 ~ r 32 T 64 ~ r 64 256 

,14 , 23,12 , 25,10 _i_ 971,8 , 421,6 i 307,4 , 107,2 49 

I f 4 -i- , -i- fi4 -1- M -i- K4 -r fi4 5Kfi 


46, 53 
52 


W W 
W 


t 4 - 7t 3 + Wt 2 - 24i + 37 


,4 1,3 ■ 55,2 19, , 1 
2 " t_ 4 8 ^ 16 

(4 . U3 , 55j2 , 19+, J_ 
2 4 816 


54 

55, 56 


w 

W\j 
7?\,7 


K 


t 3 - 5t 2 - lOt - 11 


t 3 , 11 , _ I 

1 T 4 l- g 

f3 , 11+ , I 
1 ^ 4 1 ^ 8 
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